Current large deviation function for the open asymmetric simple exclusion process 
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We consider the one dimensional asymmetric exclusion process witii particle injection and extrac- 
tion at two boundaries. Tiie model is known to exhibit four distinct phases in its stationary state. 
We analyze the current statistics at the first site in the low and high density phases. In the limit of 
infinite system size, we conjecture an exact expression for the current large deviation function. 
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Introduction. One of the main open problems in clas- 
sical statistical physics is the formulation and deriva- 
tion of simple laws that determine macroscopic quanti- 
ties in strongly interacting systems far from equilibrium. 
A broad class of nonequilibrium systems can be charac- 
terized by the presence of a macroscopic current. An 
important diagnostic tool of non-equilibrium behaviour 
is then provided by the probability distribution of cur- 
rent fluctuations. The latter is suitably represented in 
terms of its moments, which arc encoded in the current 
large deviation function (LDF). LDFs play an important 
role in the application of fluctuation theorems [l|-|3[ . Mi- 
croscopic models of interacting particles provide a use- 
ful framework for studying non-equilibrium properties in 
current-carrying classical systems and have become a ma- 
jor subject of research over the past two decades. One 
of their main uses is that their large deviation properties 
can be derived microscopically, which furnishes rigorous 
tests of underlying assumptions in phenomenological ap- 
proaches. 

The asymmetric simple exclusion process (ASEP), de- 
scribing the asymmetric diffusion of hard-core particles 
along a one-dimensional chain, is one of the best studied 
paradigms of non-equilibrium Statistical Mechanics 
The ASEP is of general interest due to its close relation 
to growth phenomena Q, as observed in recent experi- 
ments on electroconvection Q . It is also used as a model 
of molecular diffusion in zeolites 0, of biopolymers 
and sequence alignment 0, traffic flow [l^l and quan- 



FIG. 1: Dynamical rules of the ASEP. 



system evolves subject to the following rules. In the bulk 
(1 = 2,..., L~l) a particle attempts to hop one site to the 
right with rate p and one site to the left with rate q. The 
hop is executed unless the neighbouring site is occupied, 
in which case nothing happens. On the first and last sites 
these rules are modified by allowing particles to enter 
(leave) with rates a (7) at site z = 1 and with rates 5 (/?) 
at site i = L respectively, see Figure [TJ 

With every site i we associate a Boolean variable r^, 
indicating whether a particle is present (t^ = I) or not 
{Ti = 0). The state of the system at time t is then char- 
acterized by the probability distribution Pt(Ti, . . . , r^). 
The time evolution of Pt occurs according to the afore- 
mentioned rules and is subject to the master equation 
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^ = MPt. 
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tum dot chains The exact probability distribution 
for current fluctuations for the ASEP on a ring has been 
known for some time 12] ■ In the open boundary ASEP 



phenomenological [13|, approximate [Uj and numerical 
[14| treatments have been developed, but the determi- 
nation of the current LDF from first principles has been 
one of the outstanding problems in the field. Despite con- 
siderable effort, the LDF is only known in the limiting 
cases of symmetric exclusion [15[ and weak asymmetry 
[l6l |. For the infinite system the time dependence was 
obtained for total asymmetry in [l7j . 

Definition of the ASEP. At any given time t each 
site is either occupied by a particle or empty and the 



Here M = mi +771^ + mbuik is the ASEP transition ma- 
trix whose eigenvalues have non-positive real parts. The 
late time behaviour of the ASEP is dominated by the 
eigenstates of M with the largest real parts of the corre- 
sponding eigenvalues [isj . The boundary contributions 
TOi and tol describe injection (extraction) of particles 
at sites 1 and L. In the following we use a more conve- 
nient parametrization in terms of the quantities a ~ ^, 



, where 



p — g — Q! + 7± [(p — g — a + 7)^ + 4q;7] 
2^ 



(2) 



Stationary state properties of the ASEP. At late 
times the ASEP approaches a stationary state. Physi- 
cal properties then depend sensitively on the boundary 
conditions [3| • For q < p one finds four different phases 
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FIG. 2: Stationary state phase diagram for the ASEP. On the 
coexistence hne (CL) a first order phase transition occurs. 



as a function of the boundary rates as is shown in Figl2] 

Current Fluctuations. We arc interested in the prob- 
abihty distribution of the total time-integrated current 
Ql{t)^ i.e. the net number of particle jumps between 
the left boundary reservoir and site 1 in the time interval 
[0, i]. The moments of the distribution are encoded in the 
generating function (e'^'^^*^*^), where the brackets denote 
an average over all histories. In this Letter we report an 
explicit expression for the quantity 



^^(A) = lini lim - log(e 

L— >oo t— >cxo t 



AQi(t)\ 



(3) 



This characterizes the asymptotic current distribution 
which, for an ergodic system, is not expected to depend 
on the choice of initial particle configuration. 

As observed in (i^l, eqn (jS]) implies a large deviation 
property for the probability distribution P(ji,t) of the 
average current ji = Qi{t)/t at the first site. The long- 
time limiting behaviour is given by P{ji,t) ^ q-^^Ui) 
where E{ji) = ma,x\ {Xji — E{\)} is the Lcgendre trans- 
form of E{X). As a tool to compute the current LDF we 
introduce a fugacity e^ conjugate to the current on the 
first site. The boundary term mi then becomes 



mi = 



-a 
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(4) 



and E{X) is equal to the largest eigenvalue of the gener- 
alized "transition matrix" M(A). The spectrum of M( A) 
obeys a Gallavotti-Cohen symmetry [J, 0, [l^, [2^: the 
eigenvalues of Af(A') and M(A) are equal when A' and A 
are related by e^ = abcdq^^^c^^ . 

Summary of Results. Our main result is that the gen- 
erating function ^ for current fluctuations at site 1 in 



the low and high density phase and for small A and 
L — >^ cxD is of the form 



EiX) = ip-q) 



1) 



(l + a)(e^ + a)' 



(5) 



In the high density phase we obtain the same expression 
with a replaced by b. Note that the requirements that A is 
small and L — > oo explicitly break the Gallavotti-Cohen 
symmetry, as this is a duality between small and large 
negative A. We may use ([5]) to derive explicit expressions 
for the first few cumulants of the local current in terms 
of the average bulk density p = 1/(1 -I- a), 
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The first result reproduces, as expected, the bulk current 
[igj . while the second moment agrees with the diffusion 
constant in the limit g — of completely asymmetric 
diffusion [U. 

Derivation. In the following we set p ~ I without loss 
of generalit y. B ased on earlier work on the quantum XXZ 
spin chain [22| . the generalized ASEP transition matrix 
was shown to be diagonalizable using the Bethe ansatz 
in the case where the parameters satisfy ISl. |23| 



) (a/3e^ - g^/2-'=-i7<5) = 0. (7) 



Here k is an arbitrary integer in the interval |fc| < L/2. 
By considering small finite systems we find that the 
largest eigenvalue E{X) is described by one of the sets 
of Bethe equations given in [3] , which can be cast in the 
form 
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where n ~ L/2 + fc, and Yl{z) is given by 
1 



1=1 



(8) 
(9) 



9{z) + -j^9b{z) 



L 



ln(— qz) 



1=1 



Here the functions g, gb and K are given by 
(l-gz)2- 



,(z) = In 



L (1-^)2 



(10) 
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b + qz d + qz 

(12) 



K(w, z) — — ln(w) — In 



1 — qz/w 1 — q^wz 
1 — qw / z 1 — wz 



(13) 



We note that these equations are different from those 
describing the low lying excitations of the ASEP [18 1. 

The constraint (O can be satisfied for arbitrary 
a, /?, 7, 5, q and k by fixing the parameter A characterizing 
the generating function to a value among the sequences 
(SI) X'-n^ = nln(g) or (S2) A^^' ^ In (7(5g"- Va/3) , where 
n is an integer with Q < n < L. In order to infer £'(A) 
we employ the following strategy: we set A = A„ and 



then determine the ground state energies E{X^^^) of the 
corresponding generalized transition matrices. From the 
sequences of values obtained in this manner we then con- 
jecture a general expression for E{X). 

Ground State Energy for sequence (SI). Here, the 
ground state in the low density phase corresponds to a 
solution of the Bethe ansatz equations with only n roots 
(n = l,2,...) 



0( e-f"'^ 



i = 1, ■ 



(14) 



where for large L the /ij approach constant values. We 
have checked (|14p against exact diagonalization of small 
chains {L < 14) for many values of the boundary rates 
and n < 5. We conjecture that it is correct in general for 
sufficiently small n, i.e. n such that g^" > ahcdq'"^^ ~ 
'ydq^~^ /a/3. The solution is of the form of a maximal 
boundary bound state: one root lies exponentially close 
to a pole of the boundary phase shift e^'^^^ while pairs 
of the others lie on poles of the two-particle phase shift 
^K{zk,zi)_ rp]_^g ground state energy ([8]) becomes 



n-l 

3=0 



(1-9) 



a + 1 a + 



(15) 



Restoring A and p we obtain the result ([5|) . 

Ground State Energy for sequence (S2). Here the anal- 
ysis is considerably more involved. The ground state in 
the low density phase is again given by ([8]) , ([9]) , but now 
with k = L/2 — n. To keep A„ small for L ^ 1 we re- 
quire n <^ L, which corresponds to the number of Bethc 
roots being 0{L). In the following we present details for 
the case n = 1, other values can be treated analogously. 
For n = 1 there are L — 1 roots. The ground state is 
obtained by choosing 



I,^-L/2 + j, j = l,...,i-l. (16) 



The corresponding roots lie on a contour that closes as 
L — !■ oo on a point Zc on the negative real axis, see e.g. 
the plot on the left hand side of Figure |31 
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FIG. 3: Distribution of reciprocal roots 1/zj for L = 60. Left: 
a = 3.45, b = 1.5, c = -0.55, d = -0.6 and q = 0.8. Right: 
a = 1.7, b = 1.6, c = -0.55, d = -0.6 and q = 0.9. Both 
contours close on the negative real axis as L increases. 



Following 18[ we obtain an integro-differential equa- 
tion for the root density Yl{z) in the limit i — > oo, valid 
in the low and high density phases. Dropping subleading 
contributions in L"-'^ we have 



iYL{z)^g{z) + jg^,{z) + ^ 



K{w,z)Yl{w)w. (17) 



The integral from ^~ to is along the contour formed 
by the roots, and the end points are fixed by Yl{£,^) = 
±(7r — 7r/L). Equation (fT7|) may be solved by expanding 
in powers of L^^, i.e. Yl{z) = yoiz) + yi{z)/L -I- . . ., 
^ = Zc + {S + 17])/ L + . . ., which upon substitution into 
(|17p yield integro-differential equations for the functions 
2/0 and yi. Once these have been determined the corre- 

(2) 

sponding eigenvalue E{Xn ) is obtained from 
S=-^j^ e{z)Yl{z)z-^y',{zc)r]e{zc) + ..., (18) 

where we have dropped terms of 0{L~^). Here, the in- 
tegral is over the closed contour on which the roots lie. 

Assumption I: A > 0, — 1/a inside the contour. This 
regime corresponds to the case where = abed > 1, 
and is defined by assuming that —1/a lies inside the con- 
tour of integration and all other poles of g^ lie outside. 
The zeroth order term in the expansion of the counting 
function can be found as in flif , and is given by 



yo{z) 



-i In 



z 

Zc 



1- z 



(19) 



Under the above assumption, the driving term of the 
subleading integro-differential equation can be shown to 
have branch points at —1/a and at qZc ■ The branch point 
at —1/a results in branch points in yi{z) at the points 
—q^^/a, m = 0, 1, 2 . . ., and likewise for the branch point 
at qZc. As in [l8|, this suggests a functional form for 
2/1 (z) which may then be obtained explicitly. 

Finally, employing the boundary conditions for it 
is possible to show that in this regime i5 = 0, 7/ 2/g(zc) = Itt 
and that the contour closes at Zc = —bed, which agrees 
well with numerical solutions of © up to L = 200. The 
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energy can be computed from (|18p and is given by 

E^il-q)(^--^), (20) 
\a + l 1 - Zc J 

which with = —bed is fuhy consistent with ([S]), and 
coincides with it when we restore p and A using e^ = abed. 

Assumption II: A < 0, gh analytic inside the contour. 
A numerical analysis of the case e^ = abed < 1 indicates 
that the roots again lie on a contour, except for isolated 
roots on the negative real axis. The plot on the right 
hand side of Figure [3] gives an example with only one 
such isolated root zi sa —1/a. Under the assumption 
that the boundary term gi, does not have poles inside the 
contour, the leading order intcgro-difFerential equations 
may again be obtained explicitly. While the details of the 
calculation of the energy are slightly different from above, 
the final result is again (|20p with Zc = —bed, confirming 
also in this case ([5|). 

Conclusions. In this letter we have presented a con- 
jecture for the exact current LDF in the high and low 
density phases of the ASEP with open boundaries in the 
limit of infinite system size. We have presented strong 
evidence in support of our conjecture. While the den- 
sity LDF for the open ASEP has been derived from mi- 
croscopic first principles some time ago (23 |. the exact 
determination of its current LDF has been an impor- 
tant outstanding problem in non-equilibrium statistical 
physics. Both quantities are assumed to fully describe 
the experimentally accessible macroscopic behaviour of 
the ASEP . So far we have not been able to access the 
coexistence line and the maximum current phase, where 
it is necessary to scale the parameter A in a non-trivial 
way with system size [l^ . In the maximum current phase 
we have been able to analyze only the limit i — > oo for 
fixed A, in which we obtain the following result 

£;(A) = (p-q)tanh(A/4). (21) 

It would be interesting to see whether further progress 
can be made for weak asymmetry, c.f. [25|. Finally we 
note that we have obtained preliminary results on finite- 
size corrections to ([5]). 
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